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Abstract 

We present a continuum formulation of a {d + l)-dimensional directed line 
interacting with sparse potentials (i.e. (i-dimensional potentials defined only 
at discrete longitudinal locations.) An iterative solution for the partition 
function is derived. The impulsive influence of the potentials induces discon- 
tinuities in the evolution of the probability density P(x, t) of the directed line. 
The effects of these discontinuities are studied in detail for the simple case of 
a single defect. We then investigate sparse columnar potentials defined as a 
periodic array of defects in (2 + 1) dimensions, and solve exactly for P. A 
non-trivial binding/unbinding transition is found. 

PACS numbers: 03.65.-w, 05.40.+j, 68.35.Rh 
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I. INTRODUCTION 



The physics of directed lines (or directed polymers) has been the focus of much interest 
over the past two decades. This is mainly due to the very wide range of applicability of these 
simple models to important physical processes, such as wetting |T[], the motion of domain 
walls in magnets 0, and the physics of flux lines in superconductors [[J. The interest has 
been intensified by attempting to understand the effect of disordered potentials in such 
systems. In superconductors for instance, it is well known that point disorder can help to 
localise the flux lines, hence allowing true superconducting current flow; and recently, it has 
been realized that columnar defects provide an even better mechanism for pinning 0|. At 
the model level, the addition of disorder allows non-obvious connections to be made to other 
physical systems such as spin glasses [j5], non-equilibrium interface growth || and shock 
waves in fluids (in terms of the noisy Burgers equation @.) A quantitative understanding 
of the effect of such types of disorder on directed lines is still lacking. 

Aside from these important applications, the model of a single directed line interacting 
with an external potential is of substantial interest in its own right. There has been a 
great deal of work on the purely theoretical front in trying to understand these systems, 
with approaches ranging from lattice RSOS descriptions |L|, to powerful renormalization 
group studies M, and phenomenological scaling arguments ||. As indicated above, what is 
still lacking is a systematic way of treating strongly disordered potentials; although much 
progress may be made in (1+1) dimensions HlJlQfl. The overwhelming difficulties present in 
the analytic study of directed polymers may be countered to some degree by simplifying 
some aspects of the problem without trivializing the physics. One possibility is to simplify 
the form of the external potential. For instance, the problem of bulk disorder is of enormous 
interest, but completely intractable above (1+1) dimensions. One may then consider simpler 
scenarios (such as a very dilute limit of point defects, or columnar disorder) where controlled 
analytic calculations may provide precious insight. In fact, such is the richness of the physics 
of directed lines, that even for the simple example of a non-disordered columnar potential, 
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there exists non-trivial behaviour of the line, especially as a function of spatial dimension 



nun | 



In this paper we shall introduce a certain class of models which are best described as: 
'a directed line interacting with sparse potentials'. To clarify this, consider a directed line 
in (d + l)-dimensions. We label the transverse directions by a position vector y and the 
longitudinal direction by a scalar s. The line generally exists in the presence of a potential 
V(y, s). Our use of the term 'sparse potentials' corresponds to the following form for V: 

oo 
n=l 

where at this level of discussion the positions {r n } of the 'impulses' v n , along with the 
functional form of the impulses themselves are left free. Specific choices for these quantities 
may be made which then correspond to physically realizable systems. For instance, taking 
the {r n } to be regularly spaced leads one to consider a directed line interacting with a set of 
layer potentials, as might be found in a regular crystal. Alternatively, one may take the {r n } 
to be drawn from some distribution function, which along with taking v n (y) ~ S(y — y n ), 
corresponds to a directed line interacting with a very dilute set of point defects, as may be 
realized in an 'almost-pure' superconducting single crystal. 

Towards the end of the paper, we shall be interested in a specific sub-class of sparse 
potentials, namely periodic columnar arrays of point defects, where exact solutions are 
possible in the physically relevant case of (2 + l)-dimensions. Before such specialization, we 
concern ourselves with a more general analysis of the continuum theory of sparse potentials. 
An important point is that the sparse potentials have an impulsive action on the probability 
density of the directed line, leading to discontinuities in this function. Along with this effect, 
is the physical constraint that the discontinuity in the density must never be such as to make 
the density negative. The short-scale consistency of the theory at these impulses will be seen 
to have dramatic effects on the global properties of the probability density of the line. 

The outline of the remainder of the paper is as follows. In the next section, we pro- 
vide a general continuum formulation of the problem following the standard methods of a 



path-integral description. Section 3 is concerned with constructing an iterative solution for 
the partition function for a given set of sparse potentials, due account being taken of the 
discontinuous nature of this function. In section 4, we concentrate on the simple case of a 
single sparse potential which is taken to be very short-ranged in the transverse directions 
(i.e. approaching a delta-function.) This corresponds to a single point defect, and is exactly 
solvable. Simple as it may be, this case reveals the subtlety of the discontinuous nature of 
the probability density and allows us some intuition into the qualitative difference of effects 
between attractive defects (whose effect may be global) and repulsive defects (whose effect 
is always local.) In sections 5 and 6 we concentrate on the case of a periodic array of poten- 
tials, each of which is taken to be short-ranged - this is essentially a columnar array of point 
defects. Section 5 sets out the solution to the case of all the defects being either attrac- 
tive or repulsive, whilst section 6 is concerned with the solution to the case of alternating 
positive and negative defects. We restrict our attention to the physically relevant case of 
(2 + l)-dimensions, which by good fortune is the most analytically tractable. Comparison is 
made between our results for these models, and recent investigations of related microscopic 
models p!B|JT4| . We end the paper with our conclusions and a discussion of extensions to the 
present work. 



II. FORMULATION OF THE MODEL 

As mentioned in the Introduction, there are a number of model descriptions of directed 
lines, with the main difference being whether one chooses to work in the continuum or on 
a lattice. In this paper we shall use a continuum formulation which has the advantage that 
derived results will not be dependent upon microscopic parameters, and thus may be hoped 
to have some universal applicability. 

We consider a directed line in a (<i+l)-dimensional space, with y labelling the d transverse 
directions, and s labelling the longitudinal direction. We demand that the line begins at the 
point (0, 0) and ends at (x, t). The restricted partition function for the line, in the presence 
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of a potential V(y, s) is given by 

y(t)=x t t 

Z(*,t)= J Pyexp --/ 
y(0)=o I o 

where T represents temperature, and k is the elastic constant for the line. 



ds 



(i) 



Following standard methods [15] we may rewrite this path-integral in the form of a partial 



differential equation (PDE). Explicitly one finds 

d t Z(x,t) = vV 2 Z -V{*,t)Z , (2) 

with boundary condition Z(x, 0) = 5 d (x). We have defined a diffusion constant v = T/4k, 
and have scaled V so as to absorb a factor of 1/T. 

At this point we specialize to the system of interest, i.e. the case where V represents a 
set of sparse potentials. Explicitly we write 

oo 

V{y,s) = y £vMS{s-T n ) , (3) 

n=l 

where {r n } represent the longitudinal locations of the potentials {v n (y)}. As a fi- 
nal step we integrate (|2]) using the Green function of the diffusion equation g(~x.,t) = 
(<i'Kvt)~ d l 2 exp(— x 2 /4i4), which yields an integral equation for Z of the form 

Z(x, t) = g(x, t)-j d d x' J dt' g(x - x', t - t') ^ v n {x!)5(t' - r n ) Z(x', t') . (4) 



o 



71=1 



It appears to be a simple matter to integrate over the variable t' using the delta-functions, 
but as we shall see in the next section, this must be done with some care. 

We end this section with a discussion of the physical quantities which one can obtain from 
Z. It is important to realize that the restricted partition function itself is not a physical 
quantity. In order to be meaningful, it must be normalized. We therefore construct the 
probability density of the directed lines via 
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For future convenience we denote by Po( x , t), the probability density of a directed line in the 
absence of any external potential (which is actually equal to the diffusion equation Green 
function.) One may also define a local 'free energy' via 

/(x,t) = -Tlog(Z(x,t)) . (6) 

We use the term 'free energy' guardedly, since this quantity as defined above, is not extensive 
(in terms of the length t of the lines). It is, however, a useful measure of the energy /entropy 
balance for a given end-point value x. For instance, if one chooses V to be a columnar 
potential (which exists only for x = 0), then f(0,t) is a sensitive measure of a bound line 
(/ ~ t), as opposed to an unbound line (generally / ~ log(t).) In this paper we shall 
generally intuit the behaviour of the line from studying the probability density. 

As a final remark, we stress that although the path-integral and partial differential equa- 
tion forms for Z bear a striking resemblance to the Feynman path-integral and Schrodinger 
equation descriptions of a quantum mechanical particle respectively; one should use caution 
in applying results from quantum mechanics to the present problem - the idea that the 
problems are related by a simple Wick rotation (r = it) is over-simplified. Physically, the 
partition function as defined above is a positive-definite quantity. There is no analogous 
constraint on the complex wave function of quantum mechanics. Also, the time evolution 



of the wave function is generally postulated to be continuous ||16|| . In the above problem, 
there is no physical constraint on the continuity of the partition function as a function of 
t. The reason for this difference is the following. In quantum mechanics a discontinuity 
in the time evolution of the wave function corresponds to a temporal discontinuity in the 
probability density of the quantum mechanical particle, which is physically unreasonable. In 
contrast, the discontinuity of the partition function breaks no physical laws. This is because 
the symbol t used above corresponds to the length of the line - this is a fixed quantity for 
a given line. The meaning of the rate of change of the partition function with respect to t 
is to take two lines whose lengths are t and t + 8t respectively, and to equilibrate them in 
identical thermal (and possibly quenched disorder) environments; thereafter measuring the 
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partition function for each line. There is no reason a priori to insist on continuity of the 
partition function. 

These general distinctions between quantum mechanics (as expressed by Feynman's path 
integral) and the statistical mechanics of a directed line (as expressed by ([[]) above) will have 
important consequences in the sequel of this paper. 

III. DISCONTINUITIES AND 'GENERAL SOLUTION' 

To gain some insight into the nature of a sparse potential, let us simplify the problem 
to having a single potential v(x) located at a longitudinal position r. The integral equation 
for Z now takes the form 

t 

Z(x, t) = s(x, t) - j d d x' J dt' g(x. - x', t - t')v{x')5(t' - r) Z(x', t') . (7) 

o 

Obviously, for t < r we have the solution Z(x, t) = g(x, t) which implies P(x, t) = -Po( x , t) = 
g(x, t). In particular, Z~(x) = lim Z(x, r — e) = g{x, r). We now go on to the more difficult 

e— >0 

task of finding Z + (x) = limZ(x, r + e). 

In taking the limit t — > r care is needed when integrating over the delta-function in (^), 
but the Green function p(x — ~x!,t — if) = g(x. — yf,t — r) may safely be replaced by <5(x — x') 
leading to 

Z+(x) = g(x, t) - v(x) lim J dt'5{t' - r) Z(x, t') . (8) 

o 

Therefore Z + (x) = ^(x, r) +0(w), and since this also holds for Z~(x) to this order, we may 
write (H) as 

r+e 

Z+(x) = ( 7 (x,r)-^(x)limy' dt'5(t'-r) p(x,r) +0(t; 2 ) 



1 - d(x) + Off 





2a 



(9) 



It is important to note that, while one can expand the function g about t = r, this is not so 
for Z, which can now clearly be seen to be discontinuous at this point. 



To find the higher order terms in (^), we write Z + (x) and Z (x) in the unified form: 

t 

Z(x, t) = g(x, t) - v(x) J dt'5(t' - r) Z(x, £') . (10) 

o 

where t = r±e,e— >0. Iterating this equation gives 

1 - J dt'v5(t' - r) + /" d* 7 J dt"v 2 5(t' - r) <J(t" - r) + . . . I . (11) 

J 

For t < t the integrals give no contribution, confirming the result for Z~(x), and for t > r 
the lower limits of the integrals may be taken to minus infinity and t replaced by plus infinity, 
giving 



OO CO 



Z+(x) = #(x, r) i 1 - J/" dt'<u5(t' - r) + - j dt' J dt"v 2 5(t' - r) 5{f - r) + . . . 



— oo — OO 



(12) 

where the symmetry of the double integral has been used to rewrite it as an integral over 
all t and t', hence the factor of 1/2. Continuing this procedure one finds that 

00 1 r poo ~i n 

Z+(x)=s(x,r) ]T- - / dt'v5(t'-r)\ . (13) 



Therefore, in the limit e — > 0, 

U(x,r) exp(-v(x)), if t = r + e 
Z(x,t) = l (14) 

[p(x,r), ift = r-e 

From (p7Q| ) it is clear that, in the limit t — >• r, the coordinate x is simply a label and plays 
no significant part in the phenomenon we have just highlighted. Therefore, some further 
insight into this effect may be gained by a study of the zero- dimensional versions of these 
models. In order to do this it is useful to examine the differential equations corresponding 
to the above integral equations. The one corresponding to (^) has the form (cf. (|2|)): 

d t Z(x,t) = vV 2 Z -v(x)5(t-T) Z . (15) 

The zero- dimensional version of this equation is 



d t Z{t) 



vS{t~r) Z . 



(16) 



It is easy to check that the integral equation corresponding to ( jTEp is simply (|TI]) but with 



the x label absent. On the other hand, we can in this simple case, solve (|I6|) directly to find: 



in the limit e — > 0. We conclude that, while one may meaningfully formulate questions about 
the discontinuity of Z at t = r, the definition of the partition function is itself ill-defined 
precisely at this point, depending as it does on the definition of 0(0) = JL^ dtS(t). Of course, 
in a microscopic approach, the delta-function would be smoothed out and this ambiguity 
would be absent. However, as far as the evolution of Z is concerned, only its change at 
t = r, and not its actual value there, is relevant. So the conclusions that we draw in this 
paper will be independent of the precise form of any underlying microscopic model. 
Returning to the general form of the potential (|3]), the relation 



allows us to write a general iterative solution for the partition function. The idea is to split 
the evolution of Z into two parts; the first being concerned with the change in Z as the line 
encounters a sparse potential, the second with the evolution of Z between potentials. We 
naturally define 




(17) 




(18) 




(19) 



and also 




(20) 



Directly using eqn. ( f[8|) we have 




(21) 



9 



The evolution of the partition function between potentials is easily obtained since it is 
nothing but thermal wandering. We therefore have 

Z- +1 (x) = J d d x' g(x - x', r n+1 - r n ) Z+(x') . (22) 

Eqns. (^T|) and (|22|) are the main results of this section and constitute an iterative 
solution for the partition function, in some ways analogous to the usual transfer matrix 
solution used in discrete lattice formulations of directed walks [|TJ. Once the set of functions 
Z+ and Z~ is known, the partition function at intermediate values of the longitudinal co- 
ordinate may be found by quadrature from (f|). 

IV. SINGLE DEFECT 

In this section we solve perhaps the simplest example of a sparse potential, namely a 
short-ranged potential corresponding to a single point defect, located at longitudinal location 
s = t. For convenience we choose v(y) = — pA(y), where 

A(y) = (na 2 )- d/2 exp(-y 2 /a 2 ) . (23) 

The scale a is to be regarded as the shortest transverse scale in the problem, although we 
shall always need to keep a non-zero in order to regularize the theory. [Note, in the limit of 
a — > 0, the function A becomes a Dirac delta-function.] The parameter p simply represents 
the strength of the 'defect' - for p > the defect is attractive, whilst for p < 0, the defect 
is repulsive. We can now go on to calculate the probability densities. 

Since the only potential in the system is that due to the single defect, we clearly have 
Z~(x., r) = g(x, r) which implies P~(x, r) = Po( x ) r )- Integrating Z + (x, r), as given in (|14|). 
over x in order to find the appropriate normalization leads us to an expression for the prob- 
ability density of the line on the positive side of the defect. It is convenient when discussing 
this quantity to introduce a length / = {Avt) 1 ^ 2 ^> a, which is the effective transverse wan- 
dering of the line between the origin and the defect, and then to form appropriately scaled 
versions of p and / by defining p* = p/{-Ka 2 ) d ^ 2 and I* = I /a. One then finds that 
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(25) 



(26) 



Looking at the short-range form of this expression it is immediately clear that 

P+(0,r) / \ 
P (O,t) ~ \eP* + (l*) d -lj 

whilst for |x| — > oo we have 

^ + (x,r) f (ir } 

P (x,r) \ e p* + (l*y-lj 
Examination of these expressions for the relative discontinuity of the probability density 
reveals the following effects: 
Attractive weak defect p > and p = 0(a d ) 

In this case p* = 0(1), and since I* 3> 1, the short-range discontinuity of the probability 
density (as given by fl2"5])) is of order unity, whilst the long-range discontinuity (as given by 
(HP) is of negligible size. 

Attractive strong defect p > and p = 0(1) 

The situation is markedly different here: for small |x|, P + ^> Pq and for large |x|, 
P + Pq. The two probability densities become equal at some critical value of |x| which is 
much smaller than a. 
Repulsive defect p < 

The conclusions for weak repulsive defects are exactly as for weak attractive defects. For 
strong repulsive defects P + <^ Pq at short range and P + ~ Pq at long range. 

These results can be summarized by saying that at long range (in practice for |x| > a) 
weak attractive defects and all repulsive defects have no effect. But a strong attractive defect 
does have a global effect on the probability density of the directed line, at the longitudinal 
location r. 

We shall briefly consider the form of the probability density for t > r. To do this we 
write down an integral equation of the form (|7|), but with initial time r + . Using Z + (x, r), 
as given by (|i~4|), we find: 
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P (x,t) (l*) d + eP* 

where 



P(x,f) ^ (rr + F(x,t) (e^'-l) 



P(x,t) = ^— -^exp 



x 2 



4^ 2 (I-7) 



(28) 



and where we have defined ^(t) = r/t. 

Let us examine the consequences of this result for times significantly greater than r; 
i.e. we take 7 < 1. Firstly, we note that the effect of the defect upon the probability 
density is negligible, if the defect is repulsive or attractive and weak, since the quantity 
(l*) d dominates in both the numerator and denominator of ( f27|) . The situation is again 
more interesting when we consider a strong attractive defect. It is now the quantity e p * 
which dominates, at least for small enough x. Therefore, the healing of the distribution 
function at x = follows (P/Pq) x =o ~ F(0,t) ~ 1 + dr/2t. For transverse distance |x| 
being large (actually |x| 3> t(z//r) 1//2 ), the healing does not occur except at extremely large 
longitudinal distances; the relative difference in the probability density (compared to the 
free case) satisfying {P/Po) x =oo ~ (l*) d e~ p * . There will exist a scale L(t) at which the ratio 
of P to P is exactly unity, given by P(x, t) ~ 1. From Eq.(|28|) we find L(t) = (2diyt) 1 ^ 2 . 

So, to summarize the results for the simple situation of a single localized defect of strength 
p, we find that there is a qualitative difference between attractive and repulsive defects. In 
the former case, there exist two classes of defect - weak and strong - which are distinguished 
by their effect upon the probability density, this effect being local and global respectively. In 
the latter case (repulsive defect), we find that for any strength of defect, the effect upon the 
probability density is always local. The extreme asymmetry in effect of positive and negative 
localized defects will be seen to have interesting consequences in the next two sections in 
which we consider infinite arrays of defects. Since the effect of a repulsive defect upon the 
line is qualitatively the same for any strength of defect, we shall not distinguish between 
weak and strong repulsive defects. In the following sections we shall generally take the 
strength of the repulsive defect to be of order unity. 
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V. PERIODIC COLUMNAR POTENTIAL I 



In this section we shall consider a more complicated situation, namely an infinite periodic 
array of localized defects located on the column defined by x = 0. We henceforth restrict 
our attention to (2+l)-dimensions. This shares the attractive features of being both the 
most physically interesting CcLSG, clS well the most analytically tractable - a rare coincidence. 
Choosing r to be the longitudinal separation between the defects, and p to be their strength, 
we consider a potential of the form 

oo 

V(y,s) = -pY / ^(y)S(s-rn) , (29) 

n=l 

where we adopt the gaussian envelope form ( p3|) for the short-range function A. The range 
of the potentials is of 0(a) which we take to be the smallest transverse scale in the problem. 
In particular we have a <C / where I = (Aut) 1 / 2 is the effective transverse wandering of the 
line between defects. All results will be derived to leading order in a; the fact that the 
transverse scale of the function A(x) is of 0(a) allows us to frequently implement it as a 
Dirac delta function to get results to this order. We have chosen all the defects to have 
the same strength p which we can take to be either positive (attractive defects) or negative 
(repulsive defects). In the next section we shall study an analogous situation, but with 
alternating attractive and repulsive defects. 



Taking the general iterative solution as given in eqns. (|2lD and (^) and substituting the 
explicit form for the potential above, yields the relations 

Z- +1 (x) = j d d x> g(x - x\ r)Z+(x') , (30) 

and 

Z+(x) = Z"(x) exp(pA(x)) . (31) 

Combining these two results and implementing A as a Dirac delta function wherever possible 
we find 
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•7- . ^n-l(O) / \ 

Z m (x) = . ff(x,r) 



A(O) 



^n-l(O) 



e pA(o) - l] + | <fV (?(x - x', r) Z~_ x {y!) 



A(O) 

Repeating this procedure leads us to 



<?(x, r) [e^W - ll + / dV (?(x - x', 2r) Z„-_ 2 (x') exppA(x') (32) 



n-1 



( X ) = #( X > ™ T ) + R #( X ' (™ ~ m ) T ) Z m(°)- 
m=l 



(33) 



where i? is defined by 



R = 



e^°) - 1 
AfO) 



(34) 



It is convenient at this point to define ip n = Z n (0), along with f n = g(0,rn) = \/{iml 2 ). 
Setting x = in the above equation then gives 



n— 1 



A = fn + R fn-mA. 



(35) 



m=l 



This discrete equation may be solved exactly by making use of a generating function. The 
details of the calculation are relegated to Appendix A. The resulting form for ip n depends 
on whether defects are attractive or repulsive. Thus we consider these two cases separately. 

A. Attractive defects 



From Appendix A, the asymptotic (i.e. n ^> 1) result for ip n takes the form 



-1/R' 



R i(l _ e -i/R') 



exp 



n In 



1 - e-V R ' 



(36) 



where R' = R/nl 2 = [e p * — l]/(/*) 2 . So, for any p > 0, the partition function, evaluated 
at a defect site, grows exponentially; but with a rate which vanishes exponentially fast for 
small p. In order to obtain the physically meaningful probability density, we must normalize 
the partition function. Generally, the normalization is defined as N(t) = J d d xZ(-x,t). In 
particular, we define 
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N n = lim N(tti - e) . 



(37) 



Then by integrating Eq. (|33|) over the transverse space we obtain 



n-l 



N n = 1 + R ^ 



(38) 



m=l 



which on evaluating the sum gives, for large n, 



N n ~ {I /B!) exp 



n In 



1 - e-V RI 



(39) 



Dividing the partition function (|36|) by this normalization (|3^) yields the asymptotic form 
for the probability density at the n th defect: 



P B (0) S gmP(0,™-e)~ . 



(40) 



From this result we see that for any p > 0, the probability density on the column (actually 
on a defect site) attains a non-zero, constant value as n — > oo. This indicates that the line 
is always bound to the array of defects, regardless of how weakly attractive they are, or how 
widely separated. 

It is interesting to calculate the probability density on the column, but in between the 
defects. This may be directly evaluated by making use of eqn. @. Setting t = r(n + 6) 
with 6 G (0, 1] we find the asymptotic result 



P(0,t) 



R! e- 1 ^' 



R(l-e 



Z—WYeFe{-Hl-e-^' 



(41) 



where 



Fe{p) 



-u8 



du 



l-e~ u ) 



(42) 



On setting 9 — 1, the above expression reduces to the asymptotic result for defect sites, as 
given by (g[J). 

In the limit of very strong attractive defects, the probability density along the column 
has the form P(0, (n + 9)r) = (Onl 2 )^ 1 for < 6 < 1. The density takes its largest value 
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on the positive side of the defect, and then decays as 1/9 until the next defect is reached 
(0 = 1.) 

We also note that in the limit of vanishing defect strength P(0, t) reduces to the same 
form for both defect sites, and positions in between. Explicitly one has 



P(0,t) ~ e 



jp , p* < 1 , t > t . 



(43) 



In the sense that the array binds the line for any p > 0, we may say that it acts in 
precisely the same way as a constant energy column (CEC) which is attractive fi |TT| , |T2|j . 111 
that case, one has a potential of the form V(y, s) = — pA(y). The asymptotic form of the 
probability density on the column, for vanishingly small p, takes the form 



P(0,t) ~e- irp /a 2 p 



(44) 



where p = p/(47rz/). Comparing these results, we see that concerning the dominant expo- 
nential behaviour, there is an effective renormalization of the defect strength, such that it 
appears as the strength of a CEC. The precise form of this renormalization is p = p/r which 
is an intuitively appealing result. 

B. Repulsive defects 



The solution to this discrete equation for the partition function is outlined in Appendix 
A, along with details of the evaluation of the normalization. On dividing ip n by N n , we 
obtain the probability density at a defect site. We find it to have the asymptotic form 

-l 

-1 I n-rrl z 

P n (0) = \mxl 2 {B!) 2 \ln{ne 1/lR ' 1 



mrl 2 

nl 2 (R') 2 n(hi(n)) 2 



1 < n < e 1 /^' 1 

(45) 



n 



So in the deep asymptotic regime, the probability density at a defect location decays as 
Pn{0) ~ [n(ln(n)) 2 ] _1 . Again, it is interesting to compare this result to that obtained for 
the case of a CEC (this time with repulsive energy.) Following the methods of ref. one 
may ascertain that for the CEC, P(0,t) ~ [tln(t)] -1 in the asymptotic regime. It therefore 
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appears as if the defects repel the line more effectively than a CEC; which is counter-intuitive. 
The situation may be clarified by calculating the probability density in between the defects, 
i.e. taking t = r(n + 9) with n 3> 1. This may be done by making use of with the result 
that 



P((M)~~ 



1+0 -— +0 



(46) 



\n{t) J V#( ln W) 2 , 

This is the asymptotic behaviour of a free line. Therefore the repulsive array has no qual- 
itative effect upon the probability density except right at the defect positions. The decay 
of -P(0, t) for the CEC is marginally faster than a free line, but marginally slower than that 
of a line constrained to pass through a defect - this is an intuitively acceptable result. In 
contradistinction to the case of attractive defects, there is no effective renormalization of 
the defects into a CEC when they are repulsive. 



VI. PERIODIC COLUMNAR POTENTIAL II 

In the last section we have seen that there exists a great difference between an infinite 
array of attractive defects, and an infinite array of repulsive defects. In the former case, 
the line is bound to the array; and for small values of the potential energy, the array acts 
precisely in the same manner as a CEC. In the latter case, the line is oblivious to the column 
on which the array is defined, except directly at defect sites. In that case, the probability 
density is marginally reduced. There is no relation between the repulsive array, and a 
repulsive CEC. One may ask how the line acts when the array consists of both attractive 
and repulsive defects. In fact our initial motivation was to study the case of a random 
admixture of such defects along the column. However, this simple example of a disordered 
potential is extremely difficult to analyse. A simpler task is to arrange the attractive and 
repulsive defects in an alternating pattern along the column. The physics of such a system 
has attracted much interest recently |TB|,^4|, albeit in a microscopic formulation in terms 



of the RSOS model. Our presentation in this section will be rather brief as most of the 
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calculation may be constructed using the methods of the past two sections. Also, we shall 
content ourselves with only examining the gross features of this system; namely the location 
of the binding/ unbinding transition, and the qualitative behaviour of the critical properties 
of the line and the bound phase. 

Denoting the strength of the attractive defects by p > 0, and that of the repulsive defects 
by — cr < 0, we consider the set of sparse potentials 

oo oo 

V(y, s) = -pJ2 A(y)5(s - 2nr) + a £ A(y)6(s - (2n - l)r) . (47) 

n=l n=l 

Following a similar procedure to that used in the previous section we may derive a closed 
equation for the partition function at a defect site. In this case there is the minor compli- 
cation of having two types of defect, which may be easily accommodated in the following 
way. We denote the partition function at attractive and repulsive defect sites by ip!^ and ip° 
respectively: ip^ = Z^O) and ipn = ^2n-i(0)- To derive the iterative equation for these 
quantities, one uses the fundamental relations fl2l| ) and fl22|) and follows a similar procedure 
to that described in the previous section. This leads to the following simultaneous equations: 

Ti—l n 



i>n — fln + -R hn-2m1pm ~ S ^ /27i-2m+lVV 

m=l m=l 

n— 1 n— 1 

V>n = finr-i + R J2 hn-2m-l^m ~ S hn-2m^Z , (48) 
m=l m=l 

where R is given by (|34"D and 

S = ^-rr^ — • (49) 
A(0) v ; 

We shall relegate the explicit solution of these equations to Appendix B. The main result 
to emerge from this solution is the shifting of the binding/unbinding transition to a critical 
line in the (p, cr) plane which is illustrated in fig.l. The precise equation for this line is given 
in (|B12|), but the general structure takes the form 



*={"• itP<<P ' (50) 



where p c = ira 2 In 2. It is interesting to note that for p > p c the line is always bound, 
regardless of the strength of the repulsive sites. This result is in qualitative agreement with 



the recent work from the Fribourg group [ |T3| , |T4|| . This physics may have been intuitively 
expected following the analysis of the single defect, where the vast difference of effect between 
attractive and repulsive defects was examined in detail. In the terminology of section 4, we 
may say that strong attractive defects will always bind the line, whereas weak attractive 
defects require a critical strength in order to do this. There is no analogy of 'weak' and 
'strong' for repulsive defects. 

Some further details are examined in Appendix B which we shall summarize here. Firstly 
one may examine the behaviour of the line at criticality. In this case one finds that the 
probability density of the line at a defect site (attractive or repulsive) follows the asymptotic 
behaviour of a free line, namely P(0, nr) ~ 1/n. One may also study the bound state in 
which case one finds that the probability density saturates to a constant at both attractive 
and repulsive defect sites, although the density is a factor of (1/Z*) 2 smaller at the repulsive 
defect sites. One then has the picture that the although the line is bound, it really binds only 
to the attractive defects, and has a vanishingly small probability density at the repulsive 
defect sites. Having gleaned the main qualitative aspects of the alternating column, we shall 
end this section here and proceed to presenting our conclusions. 



VII. CONCLUSIONS 

In this paper we have examined a class of models described as 'a directed line in the pres- 
ence of sparse potentials'; with the understanding that a sparse potential is a <i-dimensional 
potential defined at a single longitudinal location. In section 3 we obtained a general iterative 
solution for the partition function which consisted of two pieces: free propagation between 
potentials, and a discontinuity when passing through a potential. In section 4 we considered 
in some detail a single short-ranged sparse potential which corresponds to a single defect. 
The probability density of the line was evaluated for both attractive and repulsive defects. 
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In the former case, the density was seen to undergo a global discontinuous change for strong 
defects, and the healing of the density (i.e. the relaxation to the density of a free line) was 
found to be incomplete for arbitrarily large longitudinal distances above the defect. The 
latter case of a repulsive defect was completely different in that for any strength of defect, 
the density of the line is undisturbed except within a small region about the defect. [As 
a brief aside we may relate this extreme asymmetry to the behaviour of a non-equilibrium 
interface evolving under the Kardar-Parisi-Zhang (KPZ) equation ||. Under the mapping 
between directed lines and the KPZ equation, an attractive (repulsive) defect corresponds 
to an upward (downward) force, acting for a short duration upon the surface. It has been 
previously observed |12]JT8[1 that in the strong-coupling regime of the KPZ equation, an up- 
wards force of sufficient strength may seed a large disturbance in the interface which then 
becomes effectively frozen. Alternatively, a downward force of arbitrarily large strength 
plays no role, since any disturbance it causes is quickly eradicated by the strong upward 
action of the KPZ non-linearity. The behaviour of the directed line under the influence of 
a single defect is seen to exhibit an analogous effect. The possibility of gaining intuition 
concerning the strong-coupling behaviour of the KPZ equation is a prime example of the 
usefulness of studying non-trivial, yet tractable, directed line models of the type considered 
in this paper.] 

In section 5 we studied a periodic array of defects arranged on a column, exclusively in 
(2+l)-dimensions, which is of most interest. For attractive defects, the line was found to 
be asymptotically bound, and indeed, the saturated from of the probability density along 
the column (for vanishingly weak potentials) was found to correspond to that obtained 
previously for a constant energy column (CEC), indicating that the line samples the defects 
in such a way as to renormalize their effect to be that of a CEC. For an array consisting 
of repulsive defects, we found that the density along the column is qualitatively unchanged 
from that of a free line, i.e. P ~ This result is logarithmically modified at the defect 
positions, having the form P ~ l/{t ln 2 (t)). There is no relation of these results to a repulsive 
CEC where one has P ~ l/(tln(t)), indicating that the sampling of the defects by the line 
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does not have a renormalizing, or smoothing effect. [The fact that the line may be bound 
by an array of attractive defects has a novel implication for the KPZ equation; namely, 
that a sequence of discrete upward impulses is sufficient to move the interface with non-zero 



velocity, similar to the effect of pushing with a constant force [P2"|.1 In the last section we 
examined an array consisting of alternating attractive (with strength p > 0) and repulsive 
(strength —o < 0) defects. It was found that the binding/ unbinding transition is shifted to 
a location in the (a, p) phase plane defined by the condition a = p/(l — p/p c )- This result 
is interesting as it indicates that for attractive defects of strength greater than the critical 
strength p c = 7ra 2 ln(2), the line will always be bound, regardless of the strength of the 
repulsive defects. This latter result is in accord with recent calculations on an equivalent 



microscopic RSOS model [T^JTJ]. The asymptotic behaviour of the line at criticality was 
found to be simply that of a free line (P(0,nr) ~ l/n). In the bound phase the line was 
found to be essentially bound to the attractive defects, the density at the repulsive defects 
being smaller by a factor of (1/Z*) 2 . 

We feel that the introduction of sparse potentials introduces some simplifying features 
into the study of directed lines. In the current paper we have examined probably the 
simplest form for these potentials; namely single defects, and periodic arrays of defects; 
although even for these simple periodic arrays there are many more features that may be 
examined, such as the spatial variation of the probability density away from the column, 
and the behaviour of the system in dimensions other than (2+1). One may also retain the 
simplifying nature of a periodic array of sparse potentials, but relax the condition used in 
this paper that the potentials are short-ranged in the transverse directions. For instance it 
would be of interest to study potentials which were periodic in the transverse dimensions, 
as one may then make contact to systems in which a directed line is interacting with a set 
of crystal layer potentials. Part of our motivation for examining sparse potentials was to 
see if analytic progress is possible for simple types of disorder - such as randomly placed 
defects, or regularly placed defects with random energy. We certainly regard such analyses 
as worthwhile and possible projects for the future. 
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APPENDIX A: 



In this appendix we outline the solution of the recurrence relation for the partition 
function evaluated at the defect sites, ip n . This is given by fl3"o*|): 

n-1 

4>n = fn + R fn-mi>m ■ ( A1 ) 
m=l 

The solution is most easily obtained by introducing the generating function 

CO 

$(z) s Y, * n 4>n , (A2) 

n=l 

along with a similar function f(z) defined in terms of {f n }- Summing ( |A"1| ) over n with the 
appropriate weight then gives 

**> - wMm • (A3) 

Inverting the relation ( |A^) using the calculus of residues then yields the solution 

^ = h 1%^ Z) = 2^r[^ [1 - RJ(zj] ' (M) 

where the contour C\ is a closed circle of radius 5; this circle being chosen so that no other 
singularities are enclosed bar the n th order pole at the origin. 
The function f(z) is easily evaluated to be 

oo ~n 

f(z) = (vr/ 2 )- 1 £ - = -(vr/ 2 )" 1 ln(l - z) , (A5) 

where the sum is guaranteed to converge since \z\ < S <S 1. It is the ease with which this 
sum may be evaluated in (2 + l)-dimensions which makes this case the most analytically 
tractable. Defining R' = R/(nl 2 ), we have the explicit form for ip n as 

^ n = 2mR f ^ [l + R'\n(l-z)] ' (A6) 

Examination of the integrand reveals that there exist two singularities in the complex plane 
apart from the pole at the origin. These are a branch point at z = 1 along with a simple 
pole at 
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z p = l- exp(-l/i?') = 1 - exp j--^- j . (A7) 

By cutting the contour C% on the negative real axis, we may wrap it around the rest of the 
complex plane as illustrated in fig. Al. In this way, we have 

+ / + residue at z p = . 

>C\ J cut 

We have thus replaced the essentially perturbative expression (|A6|) , by an expression which 
enables us to extract the strong-coupling behaviour, if it exists (this actually depends on 
the existence of a pole at radius 5 < \z\ < 1.) 

Whether the residue from the pole at z p dominates over the contribution from the branch 
cut, depends on the value of \z p \. It is clear from the form of the integrand that the integral 
will be dominated (for large n) by the pole if it lies within the unit circle. From ( |A7|) it 
follows that z p < 1 — exp(/*) 2 if the defects are repulsive and < z p < 1 if the defects are 
attractive. Thus for attractive defects the pole dominates and a calculation of the residue 
leads directly to fl36|). 

In the case of repulsive defects, the pole is seen to lie on the negative real axis far outside 
the circle \z\ = 1. Thus the residue at z p yields only an exponentially decaying contribution 
and is sub-dominant to the contribution from the cut. We therefore have § c ~ — / cut , which 
has the explicit form 

^ n ~ I (1 + x) n + l {(TTR'y + (l + R'\nx) 2 } ' ^ 
Referring to fl38|) we may express the normalization in integral form. Explicitly one finds 



. dx J (1 + x) n - 1 1 J 1 1 

" i T \ (1 + J j(7ri?') 2 + (l + ^lna;) 2 J ' ( } 

We shall briefly describe the asymptotic evaluation of ( |A8| ) . The integral form of the 
normalization (as well as the similar integrals which appear in the evaluation of P(0, t) for 
t 7^ nr) may be done in analogous fashion. So referring to (|A8|) , as a first step we scale x 
by n, and use the relation exp(p) = lim n ^ 00 (l +p/n) n . We then have 
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^ ~^wj dx e ~ x \** + Mx/0)] 2 I ' (A10) 

where (3 = ne~ l ^ R ' > ne^ 1 *^ ^> 1. We now split the integration range into three regions 
and estimate the order of magnitude of the integral in each region. Region i) is defined by 
< x < 1/(3, and retaining only dominant terms for small x, we find (up to prefactors) 
ip$ ~ (/31n 2 (/3)) _1 . Region ii) is defined by 1/(3 < x < (3. In this region we may drop 
ln(x) in comparison to ln(/3) which gives ~ (ln 2 (/5)) _1 . Region iii) is defined by x > (3 
in which case i/)™ ~ e~@. So clearly the contribution from region ii) dominates for large 
n. In a similar way, one may establish that the normalization has the asymptotic form of 
N n ~ 0(1) + 0(l/ln(n)). Putting these results together gives the form of the probability 
density shown in fl45|). 

APPENDIX B: 

In this appendix we outline the solution of the simultaneous iterative equations fl48|) . As 
before, it is convenient to use generating functions. Thus we define the functions 



oo 



4> E (Z) EE E ^€ , *°(Z) = E ^V? • (Bl) 
n=l n=l 

We also define 

oo oo 

f E ( Z ) EE E ^ n /2n , /°(*) EE E ^hn-X • (B2) 
n=l n=l 

Summing the iterative equations over the appropriate weight and using the above definitions 
yields the algebraic equations 

i, E {z) = f E {z) + Rf E {z)i, E {z) - Sf°{z)i,°{z) (B3) 
$°(z) = f°(z) + Rf°{z)i, E {z) - Sf E {z)i,°{z) , (B4) 

which may be readily solved, yielding the solutions 

, E(z) = f E (z) + S[f E {zf-fO{zf] 
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and 



f°(z) 



(B6) 



1 - (R - S)p{z) - RS[p{zY - p{zf} ■ 



From the particular form of /„, we also have 




1 



ln(l - z 2 ) , 



(B7) 



2tt/ 2 



and 




1 




(B8) 



2tt/ 2 



Given the definition of the generating functions, we may retrieve the original partition 
functions using 



where as before C\ is a circle enclosing the origin of small enough radius such that it encloses 
no singularities other than the pole at the origin. 

We now examine the singularity structure in the complex z plane. Firstly we note that 
there now exist two branch points, at z — ±1 which we connect to infinity with cuts along 
the real axis as shown in fig. Bl. Also, any pole which may exist within the unit circle will 
have a twin reflected through the origin due to fact that the denominators of the generating 
functions are even functions of z. As before we replace the perturbative expression by the 
non-perturbative one by cutting the contour C\ and extending around the singularities in 
the complex plane. This is illustrated in fig. Bl and leads us to the expression 



The existence of a bound state will arise only from there being a pole within the unit 
circle. We thus examine the zeroes of the denominator of the generating functions - i.e. the 
positions of the poles are solutions of 




(B9) 





(BIO) 
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where we have defined R! = R/(nl 2 ), as before, and similarly S' = S/(tt1 2 ). We analyse this 
equation by first noting that if p grows with a more slowly than a 2 , in particular if p = 0(1), 
then there is no solution for any a. On the other hand, if p grows like a 2 , or faster, then so 
do R and S and the final term in (|B10|) plays only a sub-dominant role. Therefore the poles 



are situated at ±z p , with 



Zp 



1 — exp 



1/2 

(Bll) 



(R'-S') / 

As we vary the parameters p and a, the poles exist within the unit circle only when R' > S'. 
Therefore there is a critical line in the (p, a) plane which separates the region where poles 
exist from the region where they do not, which has the equation R' = S'. A more precise 
equation for this line may be obtained by including the final term in (|B10|) in the analysis 
by substituting z p = 1 — e into the equation and solving it in the limit e — > 0. The resulting 
condition for criticality is now 

S' = — . (B12) 

l + 21n2i?' v ; 

The right-hand-side is a monotonically increasing function of R', which reaches the value 

(/*)~ 2 (corresponding to o — > oo) when p = ira 2 In 2 = p c . The general features of the critical 

line are now easy to find: a linear regime near the origin and a logarithmic approach to p c 

from below. This behavior is summarized in floPD and illustrated in fig.l. 

In order to examine the behaviour of the line at criticality, we insert the critical condition 

( [BT2D into the expressions for the generating functions. One then finds that to leading order 

ijj E (z) = f E {z) and ip°(z) = f (z), which directly gives the partition function at the defect 

sites with no integration required; i.e. we have 

tf! = l/(2mr/ 2 ) , ^ = l/((2n - 1)tt/ 2 ) . (B13) 

The appropriate normalization factors may be derived in an analogous way to that de- 
scribed in section 5. Explicitly one defines 

iV-f = lim / d d x Z(x, 2nr - e) , (B14) 

e^O J 

27 



and 

N? = lim J d d x Z(x, (2n - l)r - e) . (B15) 



n 



These functions then satisfy the relations N® = N°, and 

n— 1 n 

= l + i2£^-SZ>£. (B16) 



m=l m=l 



At criticality these functions are asymptotically constants which along with the results for the 
partition functions (|B13| ) leads to the asymptotic form of the probability density following 



P(0,tit) ~ l/n. 

In order to examine the bound state one may simply ignore the sub-dominant contri- 
butions from the cuts, and evaluate the residues of the poles at ±z p . No explicit details of 
this calculation are given here as it may easily be reconstructed from the analogous case 
examined in the previous appendix. 
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Figure Captions 



Fig. 1: The phase diagram in the (er, p) plane. The critical line separating unbound and 



bound phases in given by (|50|) . 



Fig. Al: The singularity structure in the complex z plane for the evaluation of ip(z). Also 
illustrated is the deformation (solid line) of the original contour C\ (dashed line.) 

Fig. Bl: The singularity structure in the complex z plane for the evaluation of z/> (z) and 
ijj°(z). Also illustrated is the deformation (solid line) of the original contour C\ (dashed 
line.) 
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Figure 



1 Newman and McKane 'Directed lines in sparse potentials' 




Figure Al, Newman and McKane 'Directed lines in sparse potentials' 
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Figure Bl, Newman and McKane 'Directed lines in sparse potentials' 



